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The vortex structure in superconducting stripe states is studied according to the Bogoliubov-de 
Gennes theory on the two-dimensional Hubbard model with nearest-neighbor sites pairing inter- 
action. The vortex is trapped at the outside region of the stripe line, where the superconductivity 
is weak. The superconducting coherence length along the stripe direction becomes long. There 
are no eminent low-energy electronic states even near the vortex core. These characters resemble 
the Josephson vortex in layered superconductors under a parallel field. 
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Recently, much attention has been focused on the 
stripe state of underdoped high-T c cuprates. The stripe 
state was proposed to explain the static magnetic in- 
commensurate structure observed by means of elas- 
tic neutaon, scattering experiments on La2-sSr^Cu£L 
(LSCOpBH) and Lax.6_a.Ndo.4Sra.CuO4 (LNSCO)J# 
It is considered that doped holes are localized in the 
stripe region, which contributes to the one-dimensional 
(ID) metallic conduction,- 1 and the outside region of the 
stripe is an antiferromagnetic (AF) insulator. Angle- 
resolved photoemission (ARPES) experiments were car- 
ried and the 1D-Iik_, Fermi surface in the stripe state 
was observed._#BE_> In YBa 2 Cu 3 7 _ 5 , incommensu- 
rate fluctuations that were consistent with the above 
stripe concept were reported te-m eans of inelastic neu- 
tron scattering experiments. E_'1_3'__P 

In high-Tc superconductors, it is considered that the 
low-energy electronic state around the vortex is com- 
pletely different from that nf copypxitional superconduc- 
tors. Theoretical studiesBBBBE- 1 suggest that the 
low-energy electronic state around the vortex core ex- 
tends a significant distance due to the line node of the 
d-w&ve superconducting gap in high-T c superconductors. 
Thus, we expect the zero-energy peak in the local den- 
sity of states (LDOS) at the vortex core, instead of the 
quantized energy level of the conventional s-wave pair- 
ing case. However, in the direct observation of the vor- 
tex core by scanning tunneling microscopy (STM), there 
is no eminent low-energy state around the vortex core. 
There appears only a small unexplained shoulder or iso- 
lated^^k-ajLa higher energy within the superconducting 
gap .— — f c_r These results suggest that we must con- 
sider the vortex state, including the exotic character of 
thp eJectronic state which is unique to high-T c materi- 
als.BBl_f We study the effect of the stripe state in this 
paper. 

Although there are many theoretical approaches to the 
stripe state, we base this study on the self-consistent 
Hartree-Fock (HF) theory of the Hubbard model. It is 
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believed that t 
approximation. 



stripe concept is valid beyond the HF 
We can consider the metallic stripe 
state by using the self-consistent HF theory., if we con- 
sider the realistic Fermi surface topology.— J H It can re- 
produce "the ID Fermi surface with the gap near (•§,§)" 
[see Fig. 12(d) in rcf. |2c| , which i_.suggested by the re- 
sults of ARPES experiments.Q-'B'-J It also qualitatively 
reproduces the relationship between the incommensura- 
bility and the hole density n^, including the phase tran- 
sition between the diagonal stripe in the insulator phase 
(rth < 0.05) and the vertical stripe ix. the metallic and su- 
perconducting phases (nj, > 0.05)._ , _ , B ) Thus, HF theory 
can successfully describe the stripe structure in high-T c 
cuprates as a first approximation. Therefore, we further 
investigate various phenomena in the superconducting 
state by extending this theory. Once the density of states 
(DOS) remains at the Fermi energy (i.e., metallic state), 
we can produce the superconductivity by introducing 
the pairing interaction, at least as a phenomenological 
model J_3'_J In the superconducting state of this frame- 
work, the metallic ID stripe region becomes dominantly 
superconducting, and penetrates the outside AF region. 
In this paper, we investigate the vortex state under a 
magnetic field in this superconducting state, and study 
the structure of the order parameter and the electronic 
state around the vortex. 

We begin with the conventional Hubbard model on a 
two-dimensional square lattice, and introduce the mean 
field ni tC , = (a\ a ai, a ) at the i-site, where a is a spin index 
and i = (i x ,iy)- We assume a pairing interaction V be- 
tween nearest-neighbor (NN) sites. This type of-pairing 
interaction gives the d-wave superconductivity.!-]) Thus, 
the HF Hamiltonian under a magnetic field is given by 



ft = - J_ '>-i"~>.r".i-<* + U n i-<r a i,<T a h<r 



where a 
and 



, 7 ( a i.u) is a creation (annihilation) operator, 
e represents the NN site (e = ±x,±y). The 
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transfer integral is expressed as 

7T f Ti 

t itj = U.j exp[i— / A(r) ■ dr], (2) 

with the vector potential A{r) = x r in the sym- 
metric gauge, and the flux quantum </>o- For the NN 
pairs ijj = t. For the next-NN pairs situated 

on a diagonal position on the square lattice, tij = t' . 
For the third-NN pairs, which are situated along the NN 
bond direction, ijj = t" . To reproduce the Fermi surface 
topology of LSCO, we set t' = -O.Ut and t" = O.OStJB 
The essential results of this paper do not significantly 
depend on the choice of these parameter values. 

In terms of the eigen-energy E a and the wave functions 
Uai^i), v a (ri) at the i-site, the Bogoliubov-de Gennes 
equation is given by 

\- ( K t,i,j D i,j \ ( U a {Tj) \ _ „ / u a (r t ) \ 

(3) 

where K^j = —Uj + 5i t j(Um- a - fi), Am Hf 
V^- AijSjj+e and a is an index of the eigenstateO'EP 
The self-consistent condition for the pair potential and 
the number density is given by 





=E 


"^.Xfr,)/^), (4) 


n iiT = (ot T Oi jT > 


a 


u a (n)\ 2 f(E a ), (5) 


n i,l = (°a a u) 


-Ei 

a 


w a (r,)r(l-/(B a )). (6) 


The charge density rii 


= n»,T 


+ rii t i and the spin density 



S z ,i — \( n i,] ~ n i,i)- The superconducting order param- 
eter is decomposed into d- and s-wave components at 
each site i as 

Ad,a,i — \^x,i,cr ~t~ A—x,i,a Ay^ a (7) 

A 5)(Tii = (A« ji)<7 + A_i :4 , CT + A s>!iiCT + A_ s , <iCT )/4, (8) 
with 

n Ar l+ r l+i )/2 
As i<T = A l!i+ e i0 -exp[i— / A(r)-dr], (9) 

where A iii+ e^ T = {a i+s ^ai^} and A i)i+ e,x = 

We typically consider the case of a unit cell with 24 x 24 
sites, where two vortices are accommodated. Regarding 
the stripe structure, we assume a vertical stripe of the 
incommensurability S = g, i.e., an eight-site periodic 
spin structure. It is characterized by the ordering vector 
Q = 27r(i, |— S), where the lattice constant is unity. The 
spatially averaged hole density is set to rih = 1 — ~ | 
by tuning the chemical potential /i. By introducing the 
quasimomentum of the magnetic Bloch state, we obtain 
the wave function under the periodic boundary condi- 
tion whose region covers many unit cells. We consider 
the low-temperature case T — O.Olt, where the vortex 
structure is almost independent of T . 

First, we discuss the spatial structure of the order 



parameter for the <i-wave (s-wave) superconductivity 
Ad,i = Ad,i,-\ + A d!hl (A Sji = A SiijT + A S! i,;) and the 
charge density rij by comparing the usual no-stripe case 
and the stripe case. Figure 1 presents the vortex struc- 
ture in the no-stripe case U = and V = —It. The 
amplitude |A,^| in Fig. 1(a) is suppressed near vor- 
tices. The vortex centers are located in the middle and 
at the four corners of the figure. Around the vortex core, 
the s-wave component A s ^ is induced, as shown in Fig. 
1(b), because the fourfold symmetry of the system is lo- 
cally broken where the dominant d-wave order parameter 
varies spatially. Figure 1(c) shows that the charge den- 
sity rii is slightly suppressed in the vortex core region. 
This is a charging effect of the vortex car.e_which has 
been discussed in the s-wave pairing caserj'tif We show 
that this occurs even in the d-wave pairing case by means 
of the microscopic calculation. It is suggested that the 
^ charging effect is_celated to the Hall anomaly in high-T c 
superconductors. However, to explain the experimen- 
tal results, we need the opposite sign of the charging 
to those theoretically obtained thus far. We obtain the 
opposite sign in the stripe state, as mentioned later. 

The stripe structure without a vortex was reported in 
refs. |28|, ^9] and For finite dopings from half-filling, 
the AF structure is modulated. For 6 = g, the envelope 
of the AF moment (— \y* +lv S z ,i shows the wave of the 
eight-site period along the y- (or x-) direction. Near the 
stripe line where (— l)^+ l « S z ,i changes its sign and the w- 
shift occurs, the magnetic moment is suppressed. Along 
the stripe line, the doped carriers (holes) accumulate, 
and the d-wave superconductivity \Ad,i\ is large. Outside 
of the stripe, lA^,-) is suppressed. 

Figure 2 shows the vortex structure in the stripe state 
for U = At and V = -2t. Figure 2(a) is for |A d)i |. 
Here, we show the bond-centered stripe case, where the 
center of the stripe line is on a bond. The site with 
large \Ad t i\ is in the stripe region, and that with small 
lA^j) is outside of the stripe. While lA^,-) ~ 0.24 in 
the no-stripe case in Fig. 1(a), lA^I is suppressed to 
\Ad,i\ ~ 0.15 at its maximum. This is because part of the 
DOS at the Fermi energy is already used for the stripe 
formation, and the superconductivity occurs by using the 
remaining DOS. The vortex center tends to be attracted 
to the outside of the stripe, where |A^i| is weak. We 
consider various cases for the vortex center position as 
an initial state of the iterating calculation. However, the 
vortex center is shifted to the outside of the stripe in the 
final self-consistent results. This suggests that the stripe 
functions as a line pinning center for a vortex. This is 
reasonable because the condensation energy loss required 
to create a vortex is minimal. Figure 2(a) shows that 
\Ad,%\ suppressed near the vortex core recovers quickly 
in the direction perpendicular to the stripe. The effect 
of the vortex is limited to within about three sites from 
the vortex center. However, along the parallel direction, 
|Ad,j| around the vortex core slowly recovers by using 
the maximal length of the intervortex distance. That is, 
the superconducting coherence length is short (long) in 
the perpendicular (parallel) direction. This suggests that 
the stable vortex lattice configuration is modified from 
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Fig. 1. Vortex structure in the no-stripe case for U = and V = — It. (a) [Adi[, (b) |A S j|, (c) m. The area of 24 X 24 sites is plotted, 
where vortices are located in the middle and at the four corners of the figure. 




the conventional 60° triangular lattice. However, the 
vortex core structure, which we are studying, will not be 
significantly affected by the vortex lattice configuration. 

Figure 2(b) shows the small s-wave component |A Si i| 
induced by the spatial variation of the dominant |Ad t j|. 
There, the effect of the stripe is large compared with the 
vortex effect. In the stripe region, |A a ,i| is large, and 
|A a ,t| is slightly suppressed near the vortex. As shown 
in Fig. 2(c), n% is large outside of the stripe. There, n.j 
is enhanced at the vortex core. It is important to note 
that the sign of the charging is opposite to that of the 
no-stripe case. This sign change is related to the AF 
moment of the stripe state. Also, when the AF moment 
is inducecLaround the vortex, this type of sign change is 
reported.E-$ The structure of S z i is not strongly affected 
by the vortex in our results. 

Next, we investigate the electronic structure around 
the vortex by calculating the LDOS N(E,ri) — 
N-\(E, Ti) + Ni(E,ri) at the i-site, where N^(E,ri) = 
J2 a \u a (r)\ 2 S(E - E a ) for up-spin and Ni(E,Vi) = 
J2 a \ v a( r )\ 2 S(E + E a ) for down-spin contributions. At 
the site where S z<i > 0, N r (E,r t ) > N^E,^) for 
E < Ep . In our numerical calculation, we use the deriva- 
tive f'(E) of the Fermi function instead of 5(E). In this 



case, N(E, r) corresponds to the differential tunnel con- 
ductance of STM experiments. 

The no-stripe case of Fig. 1 is shown in Fig. 3(a). 
Far from the vortex, N(E, r{) shows the typical d-wave 
superconductor's DOS at zero field. The superconduct- 
ing gap A = 2|VAd,i| ~ 0.96Z. At the vortex core, the 
superconducting gap is smeared, and a low-energy peak 
appears at E ~ 0. As for the stripe case of Fig. 2, Fig. 
3(b) shows N(E, r^) for the site within the stripe region, 
and Fig. 3(c) for the site outside of the stripe. The 
metallic stripe region has more low-energy states than 
the outside region. Far from the vortex, the LDOS is 
reduced to that of the zero-field case. It is noted that 
there appears a small gap Ai (~ 0.2i) within the d- 
wave superconducting gap A (= 2|V r A ( j ji | ~ 0.6i), be- 
cause the s-wave component is induced in the d-wave 
superconductivity- [This is also understandable from the 
ARPES results, QEtBEj-' which report that the Fermi en- 
ergy state near (f, f) vanishes as a result of the stripe 
formation. Then, the low-energy state at the gap-node 
direction (■§,■§) of cZ-wave superconductivity is absent in 
the superconducting stripe state, and a small gap opens 
in the <i-wave superconducting gap. The structure of the 
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E/t E/t E/t 

Fig. 3. Local density of states N(E,ri) at the vortex core (solid lines) and far from the vortex (dotted lines), (a) The no-stripe case 
of Fig. 1. We plot N{E,ri) at the site (i x ,i y ) = (1, 1) for the vortex core, and at the (12, l)-site far from the vortex, (b) The stripe 
region in the stripe case of Fig. 2. We plot N(E,ri) at the (l,2)-site in the vortex core, and at the (12,2)-site far from the vortex, 
(c) The outside region of the stripe in the stripe case of Fig. 2. We plot N(E,ri) at the (l,l)-site in the vortex core, and at the 
(12, l)-site far from the vortex. 



peaks above Ao is due to the stripe structure. The vor- 
tex core state is completely different from that of the 
no-stripe case. There is no eminent low-energy state. At 
the vortex core, small peaks or shoulders appear above 
the small gap Ai. The* hap; similar structures to those 
of the STM results.BBBB 

The vortex structure in the stripe state has character- 
istics in common with the Josephson vortex in layered su- 
perconductors under a parallel magnetic field. Since the 
superconductivity appears dominantly on the metallic 
ID stripe lines, the interstripe coupling between the ID 
superconducting stripes can be considered as Joscphson- 
like coupling. Also, in the layered superconductors, the 
interlayer coupling is Josephson-like coupling, and the 
vortex is trapped in the interlayer space, where the su- 
perconductivity is weak. There, the coherence length 
around the vortex is long (short) parallel (perpendic- 
ular) to the layer. While there appears a zero energy 
state (precisely speaking, there is a small gap of the or- 
der A 2 / E-p in the s-wave pairing) in the vortex core state 
for a perpendicular field, the low-energy core state does 
not appear for a parallel field.. The vortex core state has 
the energy of the order Ao.ElJ In contrast to the layered 
superconductors, the stripe line can move thermally. We 
sometimes need pinning centers in order to obtain the 
static stripe state, such as the partial substitution of Nd 
in LNSCO. The vortex core is a candidate for the pinning 
center to trap a stripe line. 

In summary, we have studied the vortex structure 
in the superconducting stripe state, based on the 
Bogoliubov-de Gennes theory. The vortex state has a 
different structure from that of the usual no-stripe case. 
The vortex is trapped in the weak superconducting re- 
gion outside of the stripe. The eminent low-energy quasi- 
particle states do not appear even at the vortex core, 
where only small peaks or shoulders appear above a small 
gap. These characteristics resemble those of the Joseph- 
son vortex in the layered superconductors under a paral- 
lel magnetic field. Although Figs. 2 and 3 represent the 
bond-centered stripe case, we observe the same type of 



behavior even in the site-centered stripe case. 

Our calculation assumes the static stripe order. Thus, 
precisely speaking, our results can be applied only to the 
static stripe state near i-filling. Apart from ^-filling, 
we must consider the fluctuation effect of the stripe line, 
which will smear the structure of the static stripe case. 
However, we expect that the main characteristics such 
as the lack of an eminent low-energy electronic state will 
remain even in the fluctuating stripe case. Thus, the 
effect of stripes should be taken into account when we 
analyze the exotic LDOS structure of the vortex state as 
observed by means of STM experiments. 
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